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One of the fundamental goals of the PHENIX experiment is to understand the structure of 
cold nuclear matter, since this serves as the initial state for heavy-ion collisions. Knowing 
the initial state is vital for interpreting measurements from heavy-ion collisions. Moreover, 
the structure of the cold nucleus by itself is interesting since it is a test-bed for our under- 
standing of QCD. In particular there is the possibility of novel QCD effects such as gluon 
saturation at low-x in the nucleus. At RHIC we can probe the behavior of gluons at low-x 
by measuring the pair cross-section of di-hadrons from di-jets in d-|-Au collisions. Our 
results show a systematic decrease in the pair cross-section as one goes to smaller impact 
parameters of the nucleus, and also as one goes to lower Bjorken x. There is a possibility 
that these interesting effects come from gluon recombination at low x in the Au nucleus. 

1 Introduction 

Deuteron-gold collisions at RHIC provide a means to explore nuclear effects on the initial-state 
parton densities in the nucleus, which is vitally important to understanding the baseline pro- 
duction for Quark-Gluon Plasma studies in heavy-ion collisions. RHIC experiments have shown 
that single inclusive hadron yields in the forward (deuteron) rapidity direction for ^/s~^ = 200 
GeV d-l- Au collisions are suppressed relative to p+p collisions [H [11|3] . The mechanism for the 
suppression has not been firmly established. Many effects have been proposed for this suppres- 
sion, such as gluon saturation [51 [S] , initial state energy loss [SJ [7] , parton recombination [5] , 
multi-parton interactions [9] , and leading and higher-twist shadowing [TUl ITT] . 

One set of measurements that might help to distinguish between the competing models 
is forward azimuthally correlated di-hadron correlation functions, which directly probe di-jet 
production through their 2— >2 back-to-back peak at A0 — tt. This technique has been used 
extensively at RHIC and is described in detail elsewhere [T2 [131 [S] • The di-hadron results 
presented here were obtained from p+p and d+Au runs in 2008 with the PHENIX detector 
and include a new electromagnetic calorimeter, the Muon Piston Calorimeter (MPC), with an 
acceptance of 3.1 < rj < 3.8 in pseudorapidity and < cj) < 2tt. 

Di-hadron measurements can probe more precise ranges of parton x in a gold nucleus than 
do single hadron probes (e.g., i?dA)- At forward rapidities, a single hadron probe will cover a 
very broad range of x, 10"'^ < xau < 0.5, thus mixing together the shadowing, anti-shadowing, 
and even EMC effects [TUj. Azimuthally correlated di-hadron measurements also enhance the 
di-jet fraction in the event selection, since one selects only the back-to-back hadrons. 
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By performing several correlation measurements with particles at different pt and rapidities, 
one can systematically scan different x ranges with an observable that is enhanced for the 
leading-order perturbative QCD component. Probing the x dependence of the effect is an 
important test since most models predict that any effects should be stronger at smaller x. 
Particles at higher pseudorapidities are produced from smaller x, so measuring hadrons from 
more forward rapidities should probe smaller x. 



2 PHENIX MPC d+Au di-Hadron Correlations 

For this analysis, back-to-back tt'^-tt'^ or hadron-Tr*^ pairs are measured with one particle at mid- 
rapidity, and the other at forward rapidity. Back-to-back cluster-Tr" pairs are also measured 
where both are in the forward rapidity region. The clusters are reconstructed from the energy 
deposit of photons in the MPC, and are estimated to be at least 80% dominated by 7r"'s, with 
the remainder coming from single photons from decays of ?7's and from direct photons. Further 
details of the analysis are available in \]A\ . 

From the pairs wc extract 
the relative yield, JdA, of corre- 
lated back-to-back hadrons in 
d+Au collisions compared to 
p+p collisions scaled with the 
average number of binary nu- 
cleon collisions (Ncou), where 
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Figure 1: JdA vs x^^^^ for 4 different centrality bins. The 
points are taken from [14 and fit with the same parametric 
function as those used for the EPS09 set of nuclear pdf's [TB]. 
The errors are the statistical and systematic errors from |14j 
added in quadrature. 



and is explained in detail in 
[T5] . JdA is simply the ana- 
log of the usual nuclear mod- 
ification factor i?dA but for 
hadron pairs. The <JdA.pp are 
the p+p or d-f Au inelastic 
cross-sections, while cr^^p^ are 
the cross-sections for di-hadron 
pair production, and is used as 
a proxy for di-jets in PHENIX. 

In Fig. [T] we have plotted 
the values of JdA versus a;;^™^ 
for four different c?-|-Au central- 
is defined 
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and can be directly measured experimentally, a;^^ should be correlated with the Bjorken x 
that is probed in the nucleus, assuming that a normal leading order (LO) perturbative QCD 
framework applies for this data. In the case of 2— 7>2 LO processes, the variable s;^™^ is lower 
than xau by the mean fragmentation fraction, (z), of the struck parton in the Au nucleus. 
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From the plot, one can see that JdA decreases with increasing centrahty, or equivalently with 
increasing nuclear thickness. The suppression also increases as one goes to lower x^™^ in the 
nucleus probed by the deuteron. 

In Fig. |2j the JdA values for three different x^^^^ are plotted versus (Ncou), the mean 
number of binary collisions, in the four centrality classes depicted in Fig. [l] One can clearly 
see from this plot a systematic decrease of JdA with greater (A'coii), as well as with decreasing 
X. The decrease is approximately linear. 



3 Discussion 

In a leading order pQCD picture, the variable JdA is 



JdA = 



'dA 



jOdA 



(1) 



for partons a+b going to outgoing jets c+d, which then fragment to hadrons with longitudinal 
fractions Zc, Zd- In the above convolutions over the parton distribution functions (/), the 
parton-parton cross-section (T, and fragmentation functions T), most of the terms are expected 
to be roughly similar between p+p and d+Au except for the nuclear gluon parton distribution 
(pdf). 

Naively, JdA might be largely dominated 
by the modification to the nuclear gluon pdf, 
since most of the events with di-hadrons at 
forward rapidities consist of a high-x parton 
from the deuteron and a low-x gluon from 
the gold nucleus. Assuming this to be true, 
one can then associate JdA with the relative 
modification of the nuclear gluon distribu- 
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Figure 2: JdA vs {NcoW 
binary collisions, for different x 



the mean number of 
^Au^ values. 



One can then interpret Fig. [2] as a system- 
atic decrease in the gluon distribution when 
one goes to the thicker parts of the nucleus, 
perhaps due to recombination of the gluons, 
and that the recombination creates a propor- 
tional decrease in the number of gluons with 
increasing number of nucleons along a line in 
the nucleus. This is schematically illustrated in Fig. [3j This decrease is stronger at lower x, 
which one might expect since the transverse size of the gluons are larger for lower x. 

If nature is kind and this data can be interpreted in terms of a simple pQCD picture, then 
this data may provide valuable information on how gluons recombine in the nucleus as a function 
of the thickness of the nucleus and Bjorken x of the gluon. Furthermore, it may be possible 
to extract Rf^, which is extremely important for understanding the quark gluon plasma since 
it forms the main ingredient for production in heavy ion collisions. However, even in a simple 
pQCD picture there have been proposed new effects in d-l-Au collisions, such as initial state 
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energy loss [6l [7] , multi-parton interactions [9] , which may complicate the interpretation of 
this data. Further theoretical studies are needed before a robust extraction of i?^" is credibly 
achieved. 




Figure 3: Schematic illustration depicting the increasing overlapping of gluons with smaller 
impact parameter due to the increasing number of nucleons along a straight line through the 
nucleus. This increases the probability for gluon recombination effects. 
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